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I INTRODUCTION 

The impact of man's activities on the atmosphere and the climatic pro- 
cesses are cur^'c-ntly important questions. To help determine the effects of 
these perturbations, comprehensive three-dimensional numerical models have 
been and are being developed for predicting the atmospheric variables. 
Although on one hand, these models require accurate representations of the 
physical processes, on the other hand they also have severe requirements for 
computational efficiency. Maintaining this delicate balance often necessi- 
tates adapting existing techniques or developing entirely new ways of treat- 
ing the processes. 

Infrared radiative transfer is an important physical process in the 
atmosphere. Although the physics of infrared radiation is fairly well 
understood, the computation of accurate heating and cooling rates can be 
expensive; thus, approximations and simplified techniques must be used. The 
schemes employed in global circulation models are apparently not directly 
applicable to stratospheric problems. Therefore, for stratospheric numerical 
models new techniques or adaptions of existing techniques must be developed. 
This report describes the infrared radiative heating model developed for 
inclusion in the NASA Ames Research Center three-dimensional stratospheric 
numerical model. 

We discuss first the background literature on both the cooling rate 

* 

calculational schemes and infrared gas data. Then, we describe the line-by- 
line radiative transfer code which, using spectroscopic gas data, computes 
infrared transmittances for given absorber amounts, temperature, and pres- 
sure. Finally, we present the infrared cooling rate computer program devel- 
oped for the Ames three-dimensional stratospheric model. 


Several points should be emphasized here. First, the ultimate cooling 
rate calculation scheme has not yet been found. The cooling rate program 
described here represents a compromise between speed, accuracy, and develop- 
mental time. However, since we do not feel that it has the optimum time/ 
accuracy characteristics, we suggest possible improvements. Second, the 
line-by-line program gives the Ames researchers the flexibility to develop 
their own hybrid forms of cooling rate calculations tailored to a variety of 
different circumstances. New arid undocumented techniques, such as Wang's 
method, could be programmed to test the characteristics of those hybrid form 
Finally* we have includL-d neither detailed descriptions that fully document 
the computer codes nor extensive examples since we view these codes as evolv 
ing. Instead, we have concentrated on presenting the rationale behind the 
code development and pointing to future improvements. 
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1 1 BACKGROUND LITERATURE SURVEY 


This section provides background data on the cooling rate calculation 
techniques and on laboratory data for the infrared radiative properties of 
atmospheric gases. 

A. COOLING RATE CALCULATION TECHNIQUES 

Calculations of atmospheric infrared heating and cooling rates have a 
long history [see Goody (1964) for a historical review]. Before large com- 
puting machines were available, graphical or tabular methods were often used 
(Elsasser, 1942; Yamamoto, 1952; Elsasser and Culbertson, 1960). As 
increasingly fast computational techniques have been developed, more exten- 
sive and detailed numerical calculations have been reported, including even 
extremely detailed line-by-line computations (Dickenson, 1972). The primary 
difficulty in performing infrared cooling calculations has been ensuring the 
inclusion of the important physical details without making the calculations 
unnecessarily burdensome and time consuming. Thus, a variety of techniques 
has been developed. 

In computing the infrared heating/cooling rate the following choices 
among procedures are available: 

> Form of the radiative transfer equation: 

- Flux-differencing 

- Direct cooling rate calculation 

- Parameterization. 

> Frequency integration: 

- Line-by-line calculation 

- Band model 

• Random or Goody model 

• Malkmus model 
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I • Exponential model 

• Empirical fit 

- Emissivities. 

> Path integration (inhomogeneous effects): 

^ - Curtis-Godsen approximation 

- Differential formulation. 

> Angle integration: 

- Diffusivity factor (exponential integral approximation) 

I - Numerical integration. 

A more detailed discuss1.;n of the cooling rate calculation literature 
is included in Appendix A. A summary of the literature is contained in 
Table 1. 

B. LABORATORY DATA 

The published data sets for the infrared radiative properties of atmo 
spheric gases are summarized in Table 2. 
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III INFRARED TRANSMISSIONS OF ATMOSPHERIC GASES 

The "first principles" approach to infrared transmission of atmospheric 
gases is to first solve the quantum mechanics equations for the states of the 
particular molecule. One can then compute the absorption coefficients from 
the transition probabilities and energy levels. In practice, this approach 
is rarely used because most molecules are too complicated to yield a solution 
to the wave equations. What car. be done is to use either experimentally 
determined or computed line position, line strength, and half width data to 
compute the absorption coefficient of the gas. This procedure of computing 
the spectral absorption coefficient of individual lines is called the line- 
by-line calculation procedure. Only recently, with the availability of 
sufficient data on atmospheric gases and with high-speed computers, has this 
technique become practical. This section of the report describes the line-by- 
line transmission computer program, which allows the user to develop a variety 
of transmission parameterizations for specific uses. 

A. LINE-BY-LINE CALCULATIONS 

The physics of infrared radiation and molecular structure is described 
by Herzberg {1945), Goody (1954), and others. Basically, an energy transi- 
tion from State 1 to State 2 gives rise to an absorption (emission) coeffi- 
cient with a particular spectral shape depending upon the dominant physical 
mechanisms. The three most common broadening mechanisms are natural, pres- 
sure (or Lorentz), and thermal (or Doppler) broadening. 

Natural broadening is actually an effect of the fundamental uncertainty 
as to the location of the line. The absorption coefficient for a line for 
natural broadening is: 



(1) 
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2 

where Is the absorption coefficient or cross section In cm /molecule, S 
Is the line strength, Is the natural broadening half width, and Oq Is 
the wavenumber location of the line center. 


Pressure (or Urentz) broadening Is an effect of the collisions between 
molecules. This effect Is proportional to the pressure, hence the name 
pressure broadening. The absorption coefficient for a line for Lorentz 
broadening Is 


K 


L 


1 

” (" ■ + “L 


( 2 ) 


where 1^^® pressure broadened half width (proportional to the total 
pressure). 


Thennal (or Doppler) broadening 1s caused by the thermal motion of the 
molecule and Is proportional to the temperature. The absorption coefficient 
for a 11f\e for Doppler broadening is 



(3) 


where «p is the Doppler broadening halfwidth. 


In the lower troposphere the pressure broadening dominates and the 
lines are essentially Lorentzian. However, In the upper stratosphere the 
pressure Is low enough that Doppler broadening becomes important*, it dom- 
inates at the very highest altitudes. Thus, there is a region where both 
Lorentz and Doppler broadening are Important. The combination of Lorentz 
and Doppler broadening, known as the Voigt profile, can be written 


K 


V 


m. 



n*) * Kp(n - n’) dn' 


(4) 


The Voigt profile expression is thus not analytic and requires the evaluation 
of the integral In Eq. (4). 
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To evaluate the absorption coefficient of the gas at ?. particular wave- 
number it is necessary to sum the absorption coefficients at that wavenumber 
from all the lines, i.e.. 


M 4 

K = 2 K^ n 
j=l 



(5) 


where ^ is the location of the center of the j-th absorption line. Although 
the procedure is straightforward, in the case of an infrared vibrational- 
rotational band there can be as many as 100,000 lines. 


In the next section, we describe the computer program that evaluates 
Eqs. (1) through (5) for Og, HgO, and COg for atmospheric conditions. 

B. OUTLINE OF THE COMPUTER CODE 


The flow diagram for the computer code is shown in Figure 1. The main 
program, TRANS, contains the line data read from the AFCRL tape. The VTABLE 
subroutine calculates a table of Voigt profile values by solving the integral 
of Eq. (4) in three di^l'ferent ways. The transmission over a given spectral 
interval is defined as 

T^„(U.T.P> 5 ^ / exp[-KU] dn . (6) 

An 


Vi'hsre U is the absorber amount, T is the temperature and P is the pressure. 

The transmission is calculated for a Voigt profile by subroutine VOIGT and 
for a Lorentz profile by subroutine LORENTZ. 

The equivalent width, the quantity most often reported in the literature, 
is defined as 


A(U,T,P) 



[1 - exp(-KU)] dn 


The equivalent width is then plotted out by routine EQWIND. 




( 7 ) 




FIGURE 1 . aOW DIAGRAM FOR THE LINE-BY-LINE CALCULATION PROCEDURE 



















The program has been tested by comparing the results for the 15 ym COg 
band with the measurements of Burch, et al. (1962), The results are gener- 
ally good but a detailed analysis has not been performed. 
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IV ATMOSPHERIC COOLING RATE CALCULATIONS 


This chapter describ^is the creation of the computer program to com^ 
^;ute atmospheric cooling rates. The underlying equations are presented 
followed by an outline of the computer code and a discussion of the 
testing results. 


A. FORMULATION. 


The starting point for the calculation of infrared radiative transfer 
is the radiative transfer equation for a plane parallel atmosphere with no 
scattering: 


>■ dT = (8) 

where u is the cosine of the zenith angle, is the radiant intensity at 
frequency v, is the absorption coefficient, is the Black Body inten- 
sity, and z is the vertical coordinate. Integrating the equation to write 
the intensity in the positive z direction yields 


l*U.v) = B„(0)T,(z.0) 


V' ' V 


*f B„(z') •pK„T,(z.z') ^ 

Jn 


V V 


where 


T^(zi»Z 2^ = exp 


f 


-pK^(z, - Zj) 




(9) 


If we now integrate Eq. (2) over the hemisphere, the radiant flux, F, 
normal to z can be written as 
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•t) ''O 


« 2irB^(0)r (z.O) 


■/ 


+ 2W B^(z')KJ;,(2.z') dz' 


( 10 ) 


where now 




f -[ 

Jn 


-K ^p(z, - z^y 

P 


pdp 




The flux can also be expressed as 


F''(z) 

V 


*27tB^(0) • T; (2.0) 


■/ 


dT' 


+ Zirl B^( 2 ‘) ^ (z.z') dz‘ 


(VI) 


The radiative cooling rate is defined 


pc, 


p\3t, 


rad 


aFM 

3Z 


(12) 


There are generally two different ways of computing Eq. (12). The 
first is simply 


lE 

9Z 


(Z) = 




2Az. 


(V3) 
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FIGURE 2. FLOW DIAGRAM FOR THE LINE-BY-LINE CALCULATION PROCEDURE 
































C. DISCUSSION 


Initial testing revealed that a straightforward, accurate evaluation 
of Eq. (14) produced cooling rates with a computational speed of roughly 
two orders of magnitude too slow for inclusion in the Three-Dimensional 
Stratospheric Model. Thus, considerable acceleration of the code is neces- 
sary. We have tried the following modifications with the listed success or 
problem as a result. 


Modification 


Result 


Tabular values of the 
absorption 

Reduction in levels from 40 

Two-grid system— one for 
cool-to-space and one for 
exchange integral 

Variation in 7,z' spacing in 
Eq. (14) 


Approximate 50% reduction 
in time 

Severe loss of accuracy 
Marginal improvement 

Approximate factor of 2 
improvement 


As the above tabulation shows, we have been able to bring the program 
to within one order of magnitude of the required speed. However, tpaeting 
the requirements of the NASA Ames Three-Dimensional Numerical Model will 
require further testing. We discuss the possible improvements in the next 
section. 
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V RECOMMENDATIONS FOR FURTHER WORK 


As discussed previously, we have developed programs that give the NASA 
Ames researchers a great deal of flexibility in calculating infrared cooling 
rates. Unfortunately, the straightforward evaluation of the cooling rate 
equation requires substantially more tifiie than is desired by the Three- 
Dimensional Stratospheric Model. Although we have succeeded in accelerating 
the code by a considerable amount, it is not certain that this approach will 
yield the desired speed/accuracy criteria. Possibly, no scheme can give the 
desired results at present. On this pessimistic note we recommend that the 
following work be undertaken: 

> Continue attempts to accelerate the code. It may be 
possible to eliminate parts of the exchange integral 
evaluation. 

> Explore new techniques. The only technique that is 
substantially different is that of Wang (Wang and 
Lacis, 1979). 

> Compute the cooling for aU; three gases and account for 
for overlap. 

> Account for nonhomogeneous paths. 

We feel that the first recommendation is the most important, for if this 
approach is successful, the second effort would be unneces^csry. The third 
and fourth recommendations are merely a matter of implementing known tech- 
niques and require no new development. 
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Appendix 


SUMMARY OF MEETINGS ON INFRARED 
HEATING/COOLING RATE CALCULATION 
at 

NASA Ames Research Center 
Moffett Field. Ca. 

5-21 July 1978 


Participants: 

J. B. Pollack 
0, B. Toon 
J. A. G nkely 
R. W. Bergstrom 


Prepared by 
R. W. Bergstrom 


28 July 1978 
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The following is at summary of the author's notes from the set of 
meetings on Infrared Heating Rate Calculations at NASA Ames Research 
Center July 5-28, 1978. Also Included are corrections and comments 
supplied by the meeting participants. The meetings were an informal 
discussion about the techniques that are currently being used and or 
being developed to compute Infrared radiative transfer. This report Is 
designed to document the opinions expressed at the meetings. 

There were 7 meetings held over a two and one-half week period. The 
first meeting was an introductory one, the next four dealt with specific 
gaseous absorbers, the sixth with numerical methods and the last was a 
potpourri. This report will give a synopsis of each meeting and then a 
sunmary and reconmended course of action as of 28 July 1978 for developing 
a radiative transfer model for inclusion In the NASA Ames three-dimensional 
numerical model. 
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Meeting 1: July 5, 1978 

Some of the unique requirements for the radiation model for the three- 
dimensional stratospheric numerical model were discussed. They are the 
following: 

> Vertical resolution from the surface to 90 km. This makes 

Doppler effects important and possible non LTE effects. 

-3 

> Severe computational efficiency requirements 'v 10 seconds 
per vertical column. 

Other uses of infrared radiation heating model (or variations of it) 
are the following; 

> Two dimensional climate model. 

> Photochemical stratospheric model . 

> Time marching one-dimensional climate model. 

Previous Ames raidation models were described. Emphasis has been 
placed on comparison to laboratory measurements. Some of the concerns 
revolve around the problems of Doppler broadening, inhomogeneous tempera- 
ture and pressure effects, minor gaseous constituents and aerosol infrared 
effects. A decision was made to devote a session's discussion to the three 
major gases of importance (C 02 » H 2 O, 0^}, a session on minor constituents 
and one on numerical problems. 
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Meeting 2; July 7, 1978 

The topic of discussion was the 15 urn CO 2 band. The discussion 
proceeded from paper to paper in a somewhat arbitrary fashion. 

1. Rodgers and Ualshaw, 1966. 

> Goody random model with quantum mechanical values for the model 
parameters. 

> Problems brought up to JAC 

1) Cooling rates do not agree with Dickinson's line by line 
calculation. 

2) The Goody random model becomes saturated at present con- 
diti ons. 

3) In general lab data is only applicable to the troposphere 
so that the stratosphere must be modeled. Then the problem 
is how to match the two regions. Line by line calculations 
are 0|C for the stratosphere but the overl aping in the 
troposphere is difficult to account for. In the troposphere 
one should compare to laboratory equivalent widths. 

t 

2. Dickinson* 1973. 

> State-of-the-art calculation. Worked in the nonoverlapping regime, 
SLinmed up equivalent widths, computed escape-to-sapce functions for 
Volght profile. Included the hotbands, isotopic bands and ten^era- 
ture dependence. Dickerson reduced the thousand or so lines in 
each of the vibration-rotation bands of COg to 10 equivalent lines, 
these equivalent lines are listed in his paper on the thermal 
structure of the atmosphere of Venus. Also the contribution of 

NLTE effects In Dickinson, 1973 are listed for each band. 

12 16 

> Non LTE effects were 20 percent 0 85 km for C O 2 . 

> Standard for comparison. 



Kuo. 1977. 

> Developed a generalized absorption coefficient. Fit a linear 
curve to the fall off with frequency (similar to the exponential 
band model of Edwards and Menard). 

> Kuo's method appears to have general applicability. Extension 
of hss method, however, requires that the temperature depen- 
dence be evaluated for the distribution of absorption coefficients. 

Pels, unpublished (Davis conference). 

> Uses a tabulated line by line calculation (Dray son). 

Pels produces for present atmospheric conditions a matrix of 
transmissions between preassigned atmospheric layers. He also 
provides derivatives of the matrix elements for + 20° K excur- 
sions from the standard temperature profiles. His method Is 
not a matrix analogue of Newtonian cooling. An example of a 
matrix analogue to Newtonian cooling Is presented in a paper 
which appeared in JAS , 34, 1402 (1977). Finally, since Pel's 
calculations are line-by-line calculations it's not clear whether 
the transmissions he produces In the troposphere match laboratory 
measurements. 

> Problems: 

1) CO 2 amount can't change. 

21 extension to troposphere. 

Stull, Wyatt, PI ass, 1962. 

> Quasi Random band model. 

> Quantum mechanical temperature effects. 

> Used In previous Ames radiation models. 

Coakely, unpublished (See p. 37). 

> Assumes strong Lorentz lines (good up to 50 km). 

> Uses temperature dependence. 

> Accounts for hot and isotopic bands. Uses overlap 
procedure of Goody and Belton. 

> Comparison to Dickinson 10 percent low. 
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7. Ramanathan» 1976. 

> Uses Cess* band model. 

> Includes only the P and R but not the Q branch of the band. 

> Accounts for Doppler broadening following Cess (1978). 
Considers pure Lorentz then pure Doppler and uses the largest 
cooling rate. 
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A. H 2 O continuum 

1. Roberts, et a1., 1976 

> analysis of transmission data 

> includes temperature effects 

2. Coffey, 1977, atmospheric measurements 

3. Stephens, 1976, calculations of IR cooling 

4. Burch, et al., unpublished laboratory data 


B. Pure rotational band 

1. Stauffer and Walsh, 1966, 14-20 urn data 

2. Palmer, 1962 data 

»»used in previous Ames radiation models 

3. Rodgers and Walshaw, 1966 

> Not tested against lab data CJAC tested it against 
Hottel's data, result poor) 

> Should be tested against Palmer's data 

4. Jurica, 1970 

> Used Wyatt, Stull and PI ass's results to compute a 
temperature dependent emissivity 

5. Coakely, unpublished 

> Used Walshaw and Rodgers tp compute a temperature 
dependent emissivity 

6. Ramanathan, 1976 

> Used emissivities published by Staley and Jurica, 1970 


C. 6.3 ym band 

1. Benedict et al., unpublished 

.> Quantum mechanical calculations 

2. Burch et al., 1962 

> lab data 

3„ Ferriso, Ludwig and Thomson, 1966 

> lab data for higher temperatures from 1-22 ym 



Meeting 4: July 12, 1978 


This meeting discussed the 9.6 pm band of ozone. 

1. McClathey. et al., 1972 

> tabulates absorption coefficient 

> used in previous Ames models 

> claims the pressure is less effective than gas amount due 
to over! aping of lines (0^ lines are very close) 

2. Aida, 1975 

> theoretical calculation 

> covers the non-overlaping line limit to full overlap. No 
Doppler effects considered. Temperature dependence found to 
be weak. 

3. Rodgers, 1969 

> Malkmus model fit to Wa Is haw's data 

> reduces correctly to non-overlapping limit (good for small 
amount, small pressure) 

4. Ramanathan, 1976 

> uses Cess' band model again 

> makes point of neglect of 9.4 pm band. However, JAC stated! 

both Rodgers and Ramanathan's results were essentially identical. 

> accounted for Doppler effects in the same way as for CO 2 

5. Dickinson, 1973 

> Ramanathan showed agreement to 10 percent 

> standard to compare against for upper levels (Aida for lower) 


Heetinq 5: July 14. 1978 


This meeting reviewed and clarified several points then discussed 
methane, nitrous oxide and aerosols. 

1 . Ramanathan's models for COg and O 3 

For COg OAC maintained that Ramanathan did not Include enough parts 
of the band. However, his results when compared to Dickinson are 10 
percent. For ozone Ramanathan's model is OK but so is Rodgers. In 
general the model Ramanathan is using asympotates correctly to the 
weak and strong line limit but not to the strong line overlap. This 
should make his results In the troposphere suspect. However, for 
ozone hlF results are 10 percent of Alda's. There was « brief discus 
Sion on Inhomogeneous path problems. 

2. N 2 O and Methane 

> Burch et al - data 

> Green et al - 1964 data 

3. Aerosols 

> Wang, Domoto, 1974 - calculaticns 

> Sucessive orders of scattering 

> Modification of optical depth 
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Meeting 6: July 17, 1978 

This meeting was devoted to numerical techniques. 

1. Angle Integration 

> compare diffusivlty factor to more accurate Integration 

2. 1-requency Integration 

> Lowest level Is 5 spectral Intervals (CO 2 * O^, S-HgO) 

> JAC accounts for overlap by assuming a water vapor atmosphere 
then perturbs It. He uses 20 spectral intervals, 50 layers 
with an accurate integration. This requires 20,000 trans- 
mission evaluations and 1-2 secs 7600 CPU time. He felt 5 
spectral intervals (emissivltles) can be made to be 1 percent 
accurate on the flux. 

3. Path Integral 

> Can use flux diverences for cooling rate calculations If use 

differential form of the equation (I.e., an equation for flux 
differences where yBdT becomes Jb " ^-1/2^’ 

essentially a Romberg integration which is very accurate (but 
slow). Hence the number of layers necessary is somewhat 
Irrelevant. Normally in GCM's they use 12 layers and a trap- 
izoidal integration. The number of transmission evaluations 
goes as the number of layers squared. Coakely's method Is 
roughly 100 times as long. His estimate of the NQAR GCM's 
code 10 percent accuracy. 

> JAC felt that in the stratosphere it is much better to evaluate 
cooling rates directly at the center of the layer rather than 
by evaluating the net flux at the top and bottom of the layer 
and taking the difference. Round off problems seriously 
compromise the latter approach. 

> For inhomogeneous paths, JBP felt that a total differential method 

Is preferable to Curtls-Godson. If * f(T, P, U) goes as (T/T^)” 
(P/Pq)*" U^ then he defines K * [ P" T^ du]^. 



JAC provided the following rebuttal; 


Coninent on problem of inhomogeneous paths. For a homogeneous 
path the mean transmission is given by 

•4v 

and presumably It can be written as 
T « f(W) 
where 

H ' (T/Tg)" (P/Pp)" 

For an Inhomogeneous path 

"W f C=/S|S^g^(v) dU] . 

•4v 

Because of differences In line profiles between homogeneous and 
Inhomogeneous paths, it need not follow that 

T * f(W) where 

H * Lf(P/Po)"^'' (T/T„)"*'‘' dU]^ . 

Stick with simple Curtls-Godson approximation when possible. 

Problem of Doppler effects. Possibilities are the following: 

> Cess's method 

> line by line 

> Break up of Kuo's coefficient into strong line regime and 
Doppler regime 
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Meeting 7: July 21, 197 8 

The last meeting discussed several general topics. 

U NCAR GCM model 

> H 2 O - Sasamorl 

> Continuum - Roberts 

> CO 2 - Coakely 
0^ * Rodgers 

The Intent Is to obtain accuracies of 1 percent on flux and 10 percent 

on cooling. There is some question if this is achieved. 

2. Coakely *s CO 2 model 

> good up to 45 km, above that the results are a factor to two compared 
to Dickinson. 

> Doppler effects should be applied separately to each band (similar to 
Rama na than) 

> There is a possibility of fitting a linear curve to the equivalent 
widths of Dickinson. This would allow a uniform modification 
instead of each band separately. 

> In the 15 ym region CO 2 bands of nonoverlapping doppler lines seem 
to have an equivalent width which is proportional to tn(U) where W 
Is the optical pathlength of C 02 * 
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SUMMARY AND RECOMMENDATIONS FOR DEVELOPING AN 
INFRARED RADIATIVE TRANSFER MODEL (RWB) 


1. COg 

> Dickinson, 1973 standard of comparison tor upper atmosphere 

> Coakely's model corrects the deficiencies uf Ramanathan's model 

2. H 2 O Continuum 

> Roberts et al. latest model 

3. H 2 O rotational and S .6 um 

> Temperature dependent emissivltles, Coakely and Jurlca seem 
promising 

4. 03 

> Either Ramanathan or Rodgers's model fits Alda calculations 

> Ramanathan’s results fit Dickinson's 

5. Numerical methods 

> Difficulty is the exchange integral evaluation (time consuming) 

> unresolved question of how to evaluate Doppler effects 

RECOMMENDATIONS FOR INFRARED CALCULATIONS FOR 3-D MODEL 

(as of 28 July 1978) 

1. COg 

> Coakely's model 

> Doppler effects as Ramanathan 

2. HgO 

> Either Coakely's or Jurica's emissivltles 

3. O 3 

> Ramanathan's model 

4. Other possibilities 

a. Pels' method 

b. Wang and Hansen s method 

5. Problems 

a-, Doppler effects ad hoc 

b. Probably will require more computer time than optimum 
(1 ml msec/ column) 
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COAKELY'S COg MODEL 

The model for the equivalent width of the 15 ym COg band is derived 
by following the procedures used to compute the equivalent widths for the 
vibration rotation bands of rigid rotor molecules (Cess and Tiwari). When 
Tines do not overlap, they absorb independently of each other and the 
equivalent width of the band is obtained by summing the equ'Walent widths 
of the lines. When the lines overlap, the equivalent width is obtained 
by allowing for overlap as described by Goody and Belton (1967). 

The equivalent width of a strong pressure-broadened absorption line 
is given by . _ 

Aj = 2r%Y>U (A.l) 

Where U is the absorber amount; y is the Lorentz half-width and Sj is the 
line strength. If lines are non-overlapping, the equivalent width of a 
vibration-rotation band like those of COg in the 15 ym region would be 
given by 

A.- <. ^ (A.l) 

where P, Q, and R represent the three branches of the band. 

are approximately 

(Ai) 
(A.l) 


In the P and R branches the strengths of the lines 
given by -Sej'- 


while in ±he Q branch they are given by ^ 

CL ^ 

s, -- pse- 
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where 




Is the band strength; B is the rotation constant (In cm’V, for CO 2 
B » 0,39 cm'b and 0 = 1,44/T (T in “K). 


If we insert (A. 3) and (A. 4) into (A,1) and then insert the result 
into (A. 2), we find that the equivalent width for the band is given by 

, s'U ' 






c 2 




where we have assumed that all lines have the same half width* y. 


The above expression applies when all lines are present, i.e., J in 
(A. 2) is allowed any integer value. Although all lines are present in 
some of the hot bands and isotopic bands in the 15 m region, in others, 
like the fundamental, every other line is missing. To obtain the equivalent 
width when every other line is missing, we restrict the sum in (A. 2) to 
even J only. With this modification we obtain an expression which is identical 
to (A. 6) but with B replaced by 4B. Thus, when every other line is missing, 
the equivalent width is given by ^ 


We notice that in the above expressions, if 'v2(T^/T)^'^^, as it would Le for 
pressure •broadened lines having the same half width, then the temperature 
dependence of y cancels that of 6, and the equivalent width for a band of 
strong non-overlapping pressure broadened lines is independent of temperature. 
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Numerical values of .the equivanent width of the 15 ym band are cal- 
culated from (A. 7} and the band strengths and isotopic abundances listed 
in Table A.1. The strength of the band is given by 







where is the abundance of the i^ isotope; r^.^ is the "line present" 
parameter listed in the table; and is the strength of the nth band. 

The half -widths of the lines broadened by N» at 300° K is taken to be 
0.71 cm atm* (Yamamoto, ). With these values, the equivalent width 
of non-overlapping lines is given by ‘^ 60/7 f'*. 

noi ( I - (S. ) • 

_ »s»/t -?W/t . p 

( / 3 . 4 / €. + 7-04 € ) O 

O 


-1 -2 

where A is in cm and U is in gm cm . The exponentials in (A. 9) allow 
for the population of the first and secorid vibrational levels at higher 
temperatures and for stimulated emission. 


When the lines overlap each other, the equivalent width of the band is 
taken to be 

(L JL(i+ ft./c ) 

In (A.IO) C is derived by fitting the epxression to laboratory measurements. 
From the measurements obtained by Burch et al. (1962) we obtain C = 48 
cm’’. 

The temperature dependence assigned to C is that for the equivalent 
width of the band in the limit of completely overlapped lines (Cess and 
Tiwari, - ). It arises because at higher temperatures, rotational levels 

with large J become populated and the width of the band grows. Laboratory 
measurements of the equivalent width, however, are usually taken at room 
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f, 


t 


t 


r 


t 


temperatures or higher, and at such temperatures the hot bands contribute 
significantly to the absorption. The measurements from which C was derived 
were taken near 300°K, and therefore, the value derived for C contains 
contributions from the hot bands. At low temperatures the hot bands are not 
present and the value of C obtained from laboratory measurements at low 
temperatures could differ from that obtained with the temperature dependence 
assigned above. 

TABLE A-1. ABUNDANCE OF CO 2 ISOTOPES AND STRENGTHS OF BANDS IN 
15 ym REGION (ORAYSON, 1967). "Lines Present" indi- 
cates whether every other line is missing, 1, or 
present, 2. 


Isotope 

Abundance 

Lines Present 


0.987 

1 

2 

c’V« 

1.1 X 10 "^ 

1 

2 

c’W® 

3.9 X 10 "^ 

2 


7.9 X 10 *^ 

2 

c’W® 

4.4 X 10 "^ 

2 


f. 
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Strength (300"K) 


Band 

Transition 

cm"^ (atm-cm)”' 

Lines Present 

1 

ooVoi^o 

194 

1 

2 

01^-02°0 

4.27 

1 

3 

01^0-10°0 

5,0 

1 

4 

01^0-02^0 

15.0 

2 

5 

02V03^0 

1.0 

1 

6 

02Vll’o 

0.022 

1 

7 

02Vll^0 

0.12 

2 

8 

02^0-03^0 

0.14 

2 

9 

02^0-03^0 

0.85 

2 

10 

10®0-03^0 

0.004 

1 

11 

03^0-04^0 

0.0042 

2 

12 

03 V 12^0 

0.0059 

2 

13 

03^0-12^0 

0.00049 

2 

14 

03^0-1 2°0 

0.014 

1 

15 

loVii^o 

0.3 

1 



EQUIVALENT WIDTH COelS^mBAND 



FIGURE A-1. EQUIVALENT WIDTHS MEASURED BY BURCH ET AL. (1962) 
(DOTS) AND COMPUTED FROM TABLE A-1 (SOLID LINE). 
The measurements shown in the figure have pressure 
broadened path lengths equivalent to the path-to- 
space pathlength in the earth's atmosphere with 
CO 2 concentrations ranging from 1/2 to 10 times 
the present concentration of 330 ppmv. 
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AND OBTAINED FROM TABLE A-1 {DASHED CURVE). The cool-to- 
space contribution to the cooling rates is shown as dots 
and indicates the importance of the path-to-space path- 
length for the cooling of the stratosphere. 
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